What are the consequences of combining nuclear and mitochondrial data for phylogenetic analysis? Lessons from Plethodon salamanders and 13 other vertebrate clades

<p>Abstract</p> <p>Background</p> <p>The use of mitochondrial DNA data in phylogenetics is controversial, yet studies that combine mitochondrial and nuclear DNA data (mtDNA and nucDNA) to estimate phylogeny are common, especially in vertebrates. Surprisingly, the consequences of combining these data types are largely unexplored, and many fundamental questions remain unaddressed in the literature. For example, how much do trees from mtDNA and nucDNA differ? How are topological conflicts between these data types typically resolved in the combined-data tree? What determines whether a node will be resolved in favor of mtDNA or nucDNA, and are there any generalities that can be made regarding resolution of mtDNA-nucDNA conflicts in combined-data trees? Here, we address these and related questions using new and published nucDNA and mtDNA data for <it>Plethodon </it>salamanders and published data from 13 other vertebrate clades (including fish, frogs, lizards, birds, turtles, and mammals).</p> <p>Results</p> <p>We find widespread discordance between trees from mtDNA and nucDNA (30-70% of nodes disagree per clade), but this discordance is typically not strongly supported. Despite often having larger numbers of variable characters, mtDNA data do not typically dominate combined-data analyses, and combined-data trees often share more nodes with trees from nucDNA alone. There is no relationship between the proportion of nodes shared between combined-data and mtDNA trees and relative numbers of variable characters or levels of homoplasy in the mtDNA and nucDNA data sets. Congruence between trees from mtDNA and nucDNA is higher on branches that are longer and deeper in the combined-data tree, but whether a conflicting node will be resolved in favor mtDNA or nucDNA is unrelated to branch length. Conflicts that are resolved in favor of nucDNA tend to occur at deeper nodes in the combined-data tree. In contrast to these overall trends, we find that <it>Plethodon </it>have an unusually large number of strongly supported conflicts between data types, which are generally resolved in favor of mtDNA in the combined-data tree (despite the large number of nuclear loci sampled).</p> <p>Conclusions</p> <p>Overall, our results from 14 vertebrate clades show that combined-data analyses are not necessarily dominated by the more variable mtDNA data sets. However, given cases like <it>Plethodon</it>, there is also the need for routine checking of incongruence between mtDNA and nucDNA data and its impacts on combined-data analyses.</p

Digits Lost or Gained? Evidence for Pedal Evolution in the Dwarf Salamander Complex (Eurycea, Plethodontidae)

Change in digit number, particularly digit loss, has occurred repeatedly over the evolutionary history of tetrapods. Although digit loss has been documented among distantly related species of salamanders, it is relatively uncommon in this amphibian order. For example, reduction from five to four toes appears to have evolved just three times in the morphologically and ecologically diverse family Plethodontidae. Here we report a molecular phylogenetic analysis for one of these four-toed lineages – the Eurycea quadridigitata complex (dwarf salamanders) – emphasizing relationships to other species in the genus. A multilocus phylogeny reveals that dwarf salamanders are paraphyletic with respect to a complex of five-toed, paedomorphic Eurycea from the Edwards Plateau in Texas. We use this phylogeny to examine evolution of digit number within the dwarf−Edwards Plateau clade, testing contrasting hypotheses of digit loss (parallelism among dwarf salamanders) versus digit gain (re-evolution in the Edwards Plateau complex). Bayes factors analysis provides statistical support for a five-toed common ancestor at the dwarf-Edwards node, favoring, slightly, the parallelism hypothesis for digit loss. More importantly, our phylogenetic results pinpoint a rare event in the pedal evolution of plethodontid salamanders

An Estimation of Erinaceidae Phylogeny: A Combined Analysis Approach

BACKGROUND: Erinaceidae is a family of small mammals that include the spiny hedgehogs (Erinaceinae) and the silky-furred moonrats and gymnures (Galericinae). These animals are widely distributed across Eurasia and Africa, from the tundra to the tropics and the deserts to damp forests. The importance of these animals lies in the fact that they are the oldest known living placental mammals, which are well represented in the fossil record, a rarity fact given their size and vulnerability to destruction during fossilization. Although the Family has been well studied, their phylogenetic relationships remain controversial. To test previous phylogenetic hypotheses, we combined molecular and morphological data sets, including representatives of all the genera. METHODOLOGY AND PRINCIPAL FINDINGS: We included in the analyses 3,218 bp mitochondrial genes, one hundred and thirty-five morphological characters, twenty-two extant erinaceid taxa, and five outgroup taxa. Phylogenetic relationships were reconstructed using both partitioned and combined data sets. As in previous analyses, our results strongly support the monophyly of both subfamilies (Galericinae and Erinaceinae), the Hylomys group (to include Neotetracus and Neohylomys), and a sister-relationship of Atelerix and Erinaceus. As well, we verified that the extremely long branch lengths within the Galericinae are consistent with their fossil records. Not surprisingly, we found significant incongruence between the phylogenetic signals of the genes and the morphological characters, specifically in the case of Hylomys parvus, Mesechinus, and relationships between Hemiechinus and Paraechinus. CONCLUSIONS: Although we discovered new clues to understanding the evolutionary relationships within the Erinaceidae, our results nonetheless, strongly suggest that more robust analyses employing more complete taxon sampling (to include fossils) and multiple unlinked genes would greatly enhance our understanding of the Erinaceidae. Until then, we have left the nomenclature of the taxa unchanged; hence it does not yet precisely reflect their phylogenetic relationships or the depth of their genetic diversity

Molecular evolution of HoxA13 and the multiple origins of limbless morphologies in amphibians and reptiles

Developmental processes and their results, morphological characters, are inherited through transmission of genes regulating development. While there is ample evidence that cis-regulatory elements tend to be modular, with sequence segments dedicated to different roles, the situation for proteins is less clear, being particularly complex for transcription factors with multiple functions. Some motifs mediating protein-protein interactions may be exclusive to particular developmental roles, but it is also possible that motifs are mostly shared among different processes. Here we focus on HoxA13, a protein essential for limb development. We asked whether the HoxA13 amino acid sequence evolved similarly in three limbless clades: Gymnophiona, Amphisbaenia and Serpentes. We explored variation in ω (dN/dS) using a maximum-likelihood framework and HoxA13sequences from 47 species. Comparisons of evolutionary models provided low ω global values and no evidence that HoxA13 experienced relaxed selection in limbless clades. Branch-site models failed to detect evidence for positive selection acting on any site along branches of Amphisbaena and Gymnophiona, while three sites were identified in Serpentes. Examination of alignments did not reveal consistent sequence differences between limbed and limbless species. We conclude that HoxA13 has no modules exclusive to limb development, which may be explained by its involvement in multiple developmental processes

Galaxies and Cladistics

The Hubble tuning fork diagram, based on morphology and established in the 1930s, has always been the preferred scheme for classification of galaxies. However, the current large amount of multiwavelength data, most often spectra, for objects up to very high distances, asks for more sophisticated statistical approaches. Interpreting formation and evolution of galaxies as a ?transmission with modification' process, we have shown that the concepts and tools of phylogenetic systematics can be heuristically transposed to the case of galaxies. This approach, which we call ?astrocladistics', has successfully been applied on several samples. Many difficulties still remain, some of them being specific to the nature of both galaxies and their diversification processes, some others being classical in cladistics, like the pertinence of the descriptors in conveying any useful evolutionary information.Comment: Talk given at the "12th Evolutionary Biology Meeting" held in Marseille, France, Sept. 24-26, 200

Macro-Climatic Distribution Limits Show Both Niche Expansion and Niche Specialization among C4 Panicoids

Grasses are ancestrally tropical understory species whose current dominance in warm open habitats is linked to the evolution of C4 photosynthesis. C4 grasses maintain high rates of photosynthesis in warm and water stressed environments, and the syndrome is considered to induce niche shifts into these habitats while adaptation to cold ones may be compromised. Global biogeographic analyses of C4 grasses have, however, concentrated on diversity patterns, while paying little attention to distributional limits. Using phylogenetic contrast analyses, we compared macro-climatic distribution limits among ~1300 grasses from the subfamily Panicoideae, which includes 4/5 of the known photosynthetic transitions in grasses. We explored whether evolution of C4 photosynthesis correlates with niche expansions, niche changes, or stasis at subfamily level and within the two tribes Paniceae and Paspaleae. We compared the climatic extremes of growing season temperatures, aridity, and mean temperatures of the coldest months. We found support for all the known biogeographic distribution patterns of C4 species, these patterns were, however, formed both by niche expansion and niche changes. The only ubiquitous response to a change in the photosynthetic pathway within Panicoideae was a niche expansion of the C4 species into regions with higher growing season temperatures, but without a withdrawal from the inherited climate niche. Other patterns varied among the tribes, as macro-climatic niche evolution in the American tribe Paspaleae differed from the pattern supported in the globally distributed tribe Paniceae and at family level.Fil: Aagesen, Lone. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Botánica Darwinion. Academia Nacional de Ciencias Exactas, Físicas y Naturales. Instituto de Botánica Darwinion; ArgentinaFil: Biganzoli, Fernando. Universidad de Buenos Aires. Facultad de Agronomía. Departamento de Métodos Cuantitativos y Sistemas de Información; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Bena, María Julia. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Botánica Darwinion. Academia Nacional de Ciencias Exactas, Físicas y Naturales. Instituto de Botánica Darwinion; ArgentinaFil: Godoy Bürki, Ana Carolina. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Botánica Darwinion. Academia Nacional de Ciencias Exactas, Físicas y Naturales. Instituto de Botánica Darwinion; ArgentinaFil: Reinheimer, Renata. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Instituto de Agrobiotecnología del Litoral. Universidad Nacional del Litoral. Instituto de Agrobiotecnología del Litoral; ArgentinaFil: Zuloaga, Fernando Omar. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Botánica Darwinion. Academia Nacional de Ciencias Exactas, Físicas y Naturales. Instituto de Botánica Darwinion; Argentin

Latitudinal diversity gradients in New World bats: are they a consequence of niche conservatism?

The increase in species diversity from the Poles to the Equator is a major biogeographic pattern, but the mechanisms underlying it remain obscure. Our aim is to contribute to their clarification by describing the latitudinal gradients in species richness and in evolutionary age of species of New World bats, and testing if those patterns may be explained by the niche conservatism hypothesis. Maps of species ranges were used to estimate species richness in a 100 x 100 km grid. Root distances in a molecular phylogeny were used as a proxy for the age of species, and the mean root distance of the species in each cell of the grid was estimated. Generalised additive models were used to relate latitude with both species richness and mean root distance. This was done for each of the three most specious bat families and for all Chiroptera combined. Species richness increases towards the Equator in the whole of the Chiroptera and in the Phyllostomidae and Molossidae, families that radiated in the tropics, but the opposite trend is observed in the Vespertilionidae, which has a presumed temperate origin. In the whole of the Chiroptera, and in the three main families, there were more basal species in the higher latitudes, and more derived species in tropical areas. In general, our results were not consistent with the predictions of niche conservatism. Tropical niche conservatism seems to keep bat clades of tropical origin from colonizing temperate zones, as they lack adaptations to survive cold winters, such as the capacity to hibernate. However, the lower diversity of Vespertilionidae in the Neotropics is better explained by competition with a diverse pre-existing community of bats than by niche conservatism.MJRP was supported by Foundation for Science and Technology, Portugal (www.fct.pt), fellowship SFRH/BD/19620/2004 and SFRH/BPD/ 72845/2010. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We thank J. P. Granadeiro and R. Lemos for advice on statistical methods, J. A. Diniz-Filho for support with SAM software, and M. A. Dias and P. E. Cardoso for help with the distributional maps. O. R. P. Bininda-Emonds kindly gave us access to the Chiroptera species-level phylogeny. Bat data was provided by NatureServe in collaboration with Bruce Patterson, Wes Sechrest, Marcelo Tognelli, Gerardo Ceballos, the Nature Conservancy-Migratory Bird Program, Conservation International-CABS, World Wildlife Fund-US, and Environment Canada-WILDSPACE - http://www.natureserve.org/infonatura.publishe

Macroevolutionary Patterns in the Aphidini Aphids (Hemiptera: Aphididae): Diversification, Host Association, and Biogeographic Origins

, the most diverse genus in the family. We used a combined dataset of one nuclear and four mitochondrial DNA regions. A molecular dating approach, calibrated with fossil records, was used to estimate divergence times of these taxa.Most generic divergences in Aphidini occurred in the Middle Tertiary, and species-level divergences occurred between the Middle and Late Tertiary. The ancestral state of host use for Aphidini was equivocal with respect to three states: monoecy on trees, heteroecy, and monoecy on grasses. The ancestral state of Rhopalosiphina likely included both heteroecy and monoecy, whereas that of Aphidina was most likely monoecy. The divergence times of aphid lineages at the generic or subgeneric levels are close to those of their primary hosts. The species-level divergences in aphids are consistent with the diversification of the secondary hosts, as a few examples suggest. The biogeographic origin of Aphidini as a whole was equivocal, but the major lineages within Aphidina likely separated into Nearctic, Western Palearctic, and Eastern Palearctic regions.Most generic divergences in Aphidini occurred in the Middle Tertiary when primary hosts, mainly in the Rosaceae, were diverging, whereas species-level divergences were contemporaneous with diversification of the secondary hosts such as Poaceae in the Middle to Late Tertiary. Our results suggest that evolution of host alternation within Aphidini may have occurred during the Middle Tertiary (Oligocene) when the secondary hosts emerged

Distinct Origin of the Y and St Genome in Elymus Species: Evidence from the Analysis of a Large Sample of St Genome Species Using Two Nuclear Genes

Previous cytological and single copy nuclear genes data suggested the St and Y genome in the StY-genomic Elymus species originated from different donors: the St from a diploid species in Pseudoroegneria and the Y from an unknown diploid species, which are now extinct or undiscovered. However, ITS data suggested that the Y and St genome shared the same progenitor although rather few St genome species were studied. In a recent analysis of many samples of St genome species Pseudoroegneria spicata (Pursh) À. Löve suggested that one accession of P. spicata species was the most likely donor of the Y genome. The present study tested whether intraspecific variation during sampling could affect the outcome of analyses to determining the origin of Y genome in allotetraploid StY species. We also explored the evolutionary dynamics of these species.Two single copy nuclear genes, the second largest subunit of RNA polymerase II (RPB2) and the translation elongation factor G (EF-G) sequences from 58 accessions of Pseudoroegneria and Elymus species, together with those from Hordeum (H), Agropyron (P), Australopyrum (W), Lophopyrum (E(e)), Thinopyrum (E(a)), Thinopyrum (E(b)), and Dasypyrum (V) were analyzed using maximum parsimony, maximum likelihood and Bayesian methods. Sequence comparisons among all these genomes revealed that the St and Y genomes are relatively dissimilar. Extensive sequence variations have been detected not only between the sequences from St and Y genome, but also among the sequences from diploid St genome species. Phylogenetic analyses separated the Y sequences from the St sequences.Our results confirmed that St and Y genome in Elymus species have originated from different donors, and demonstrated that intraspecific variation does not affect the identification of genome origin in polyploids. Moreover, sequence data showed evidence to support the suggestion of the genome convergent evolution in allopolyploid StY genome species